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Abstract 

The 58 iVi + 58 Ni reaction at 30 MeV/nucleon has been experimentally investigated 
at the Superconducting Cyclotron of the INFN Laboratori Nazionali del Sud. In 
midperipheral collisions the production of massive fragments (4<Z<12), consistent 
with the statistical fragmentation of the projectile-like residue and the dynamical 
formation of a neck, joining projectile-like and target-like residues, has been ob- 
served. The fragments coming from these different processes differ both in charge 
distribution and isotopic composition. In particular it is shown that these mecha- 
nisms leading to fragment production act contemporarily inside the same event. 
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The production of intermediate mass fragments (IMF, 3<Z<20) is the distinguishing feature 
of intermediate energy heavy ion collisions. A possible scenario for IMF emission involves 
the development of bulk instabilities in the nuclear matter, which lead to a statistical frag- 
mentation. This is supported by several experimental results concerning both the decay of 



unique sources formed in central collisions and the disassembly of quasitarget (QT) and 
quasiprojectile (QP) in dissipative peripheral collisions [1,2]. 

On the other hand theoretical studies of the collision dynamics predict the formation of 
highly deformed nuclei in this energy regime [3] and experimental evidences have been col- 
lected on the role of the formation and decay of neck-like structures in IMF production [4,5]. 
Nevertheless it is still unclear when this reaction mechanism sets in. For medium mass sys- 
tems Boltzmann-Nordheim-Vlasov (BNV) calculations predict the onset of neck instabilities 
to occur for energies between 30 and 70 MeV/nucleon [3]. For lower energies the interaction 
time between the projectile and the target is sufficiently long for the neck to be reabsorbed 
by the QP and/or the QT, whereas for higher energies the system evolves towards a fireball 
regime. 

The study of the formation and decay of neck-like structures has a variety of theoretical 
implications. The beam energy at which a neck zone is formed and decays depends on 
the Equation of State (EOS) of the nuclear matter [6]; a soft EOS favours neck formation 
(because of a smaller orbiting) and rupture. Similarly an isospin dependence of the EOS is 
also expected to influence the process of formation and decay of an intermediate source and 
this should reflect on the neutron-to-proton ratio of the emitted fragments [7] . 

The present analysis of the Ni+Ni midperipheral collisions is aimed to investigate the inter- 
play, inside the same event, between dynamically driven neck instabilities and QP statistical 
fragmentation in IMF production. 

The experiment was performed at the INFN Laboratori Nazionali del Sud with MEDEA [8] 
and MULTICS [9] apparata. A beam of 58 Ni at 30 MeV/nucleon bombarded a 2 mg/cm 2 
thick nickel target. The angular range 3° < 9[ a b <28° was covered by the MULTICS array 
[9], consisting of 55 telescopes, each of which was composed of an Ionization Chamber (IC), 
a Silicon detector (Si) and a Csl crystal. Typical energy resolutions were 2%, 1% and 5% 
for IC, Si and Csl, respectively. The threshold for charge identification in the MULTICS 
array was about 1.5 MeV/nucleon. A good mass resolution for Z=l-6 isotopes was obtained 
above 8.5, 10.5, 14 MeV/nucleon for 4 He, 6 Li and 12 C nuclei, respectively [2]. Light charged 
particles (Z=l,2) and 7-rays were detected at 30° < 9 tab <170° by the BaF 2 ball of the 
MEDEA apparatus. The geometric acceptance of the combined array was greater than 90% 

Of 47T. 

Since the relatively small number of detected fragments and particles (Nc<5-6) makes not 
possible an accurate impact parameter selection, using the standard methods [10], we fol- 
lowed a different approach: 1) selecting only the "complete" multifragment events, i.e. events 
where at least three IMF were produced and at least 80% of the total linear momentum was 
detected; 2) defining peripheral and midperipheral collisions those for which the heaviest 
fragment (with charge at least 1/3 of that of the projectile, i.e. Z>9) travels, in the lab- 
oratory frame, with a velocity higher than 80% of that of the projectile (vp=7.6 cm/ns). 
Accordingly, since the energy thresholds make not possible the detection of the QT reaction 
products, we find that the total detected charge (Z^ot) does not differ from that of the 
projectile for more than 30% (20<Z TOT <36). Considering that not completely detected 
central events can simulate complete non-central ones, we checked the existence in central 
events of heavy fragments (Z>9) moving faster than 5.6 cm/ns, i.e. compatible with the 
QP velocity. Central collision events have been selected requiring that the heaviest fragment 
has a velocity close to the center of mass one and that the total detected charge is at least 



80% of the initial (projectile + target) value. We found that such a kind of contamination 
is completely absent. 

To get information about the impact parameters range selected by this procedure and by the 
apparatus acceptance, the experimental data have been compared with the predictions of 
classical molecular dynamics (CMD) calculations [11]. CMD events are plotted in Fig. la; the 
dot-filled area refers to the amount of "complete" multifragmentation events (with at least 
three IMF products) after experimental efficiency filtering. In the dark area of Fig. la we 
present the filtered complete events with the further condition that the heaviest fragment 
moves with a velocity higher than the 80% of v P . It becomes evident how midperipheral 
impact parameter are preferentially selected. 




Fig. 1. (a) Yields from CMD calculations: raw data (blank area), multifragmentation events after 
apparatus efficiency (dot-filled area) , with further constraints on the heaviest fragments (dark area) , 
(b-d) Angular distributions for IMF forward emitted by the quasiprojectile (all charged particle 
(b), Z=4 (c) and Z=6 (d) fragments). 

The results presented hereafter will refer only to midperipheral events with at least three 
IMF, with the aim of observing the IMF emitting sources. To this respect we present in 
Fig. 2 the distributions of the parallel component of the velocity, with respect to the beam 
direction, for different Z values. For carbon and oxygen nuclei two distinct distributions 
are evident: the first is centered at 6.5 cm/ns (the QP velocity) and the second one is at 
the center of mass velocity (3.8 cm/ns), intermediate between that of the target and of 
the projectile, because of the system symmetry. Figg.2b-c show, in events where both the 
heaviest fragment and an oxygen nucleus have velocities close to that of the projectile, the 
yield of a third fragment (dot-dashed lines) . It appears that this latter fragment has a great 
probabilty of being emitted with a velocity centered around that of the cm.. This is evidence 
that, while the QP is decaying in two or more fragments, there is another system emitting 



IMF, intermediate (IS) between the QP and QT. 




Fig. 2. v par distributions for Z=3, 6 (dot-dashed line, multiplied by a factor 10), 8 (dot-dashed 
line, multiplied by a factor 40), 12; the arrows refer to the center of mass (CM) and QP velocities. 

With the same event selection used for the fragment analysis for midperipheral collisions 
we considered the proton spectra detected by the MEDEA apparatus. Eight energy spectra, 
measured at different laboratory angles, were simultaneously fitted by the superposition of 
three Maxwellian distribution. The fit results put in evidence the existence of a fast source 
(6.8 cm/ns and slope parameter T s i ope ~5.1 MeV) consistent with the QP, an almost at rest 
source (QT system, 0.8 cm/ns, T s i op e —4.3 MeV) and an intermediate velocity source (3.8 
cm/ns) with high slope parameter (T s i op e —10.1 MeV). Moreover the QP and QT source 
show the same proton multiplicity. The intermediate velocity component has been generally 
interpreted as also due to pre-equilibrium emission in the interacting zone [12]. These results 
are compatible with those shown in Ref. [13]. 

These results show that in midperipheral collisions three different emitting sources are 
present; there are events in which the IMF can be simultaneously produced by the decay 
of QP and QT (its fragments are not seen because under energy threshold for identifica- 
tion) sources and from a neck, forming a midvelocity emission source. Thus, disentangling 
the contributions from the QP and the neck sources becomes a mandatory requirement to 
improve the understanding of the IMF production mechanism and perform comparisons of 
the IMF experimental yields with theoretical predictions. To this purpose we first study 
the process leading to the disassembly of the QP restricting the analysis to the fragments 
emitted with v par >6.5 cm/ns (see Fig.2). This constraint allows the selection of the decay 
products forward emitted in the QP decay with negligible contamination due to QT and 
midrapidity source emissions. To check if the QP reaches an equilibration stage before its 
de-excitation we studied in its reference frame the angular and energy distributions of the 
emitted fragments. Angular distributions are presented in Figg.lb-d; they are flat and then 
in agreement with the hypothesis of an isotropic emission; this is a necessary condition to 



establish a possible equilibration of the studied system. On the other hand fitting the en- 
ergy distributions of different isotopes with Maxwellian functions we get, for all the detected 
isotopes (3<A<14), similar values (within errors) for the apparent temperatures T s i ope . This 
behaviour gives indications that the condition of equilibration of the fragmenting systems 
is satisfied. In Figg.3a-e the energy distributions of different isotopes are presented. Each 
panel (from a to d) refers to two different isotopes of the same element to put in evidence 
the similarity of shapes and slopes; in the (e) panel the Maxwellian fit is superimposed to 
the experimental data showing the accuracy of the obtained results. 

To gain more insight in the equilibration of the QP we compared the experimental charge dis- 
tribution with the microcanonical Statistical Multifragmentation Model (SMM) predictions 
[15], performed for a Ni nucleus at one third of the normal density. The events generated 
by SMM at different input excitation energies were filtered by the experimental acceptance. 
The experimental charge distribution is quite well reproduced assuming an excitation energy 
of 4 MeV/nucleon (Fig.3f). 




Fig. 3. (a-e) Energy distributions for different isotopes (a - 12 C (full line), 13 C (dashed line), b - 
10 B (full), n B (dashed), c - 7 Be (full), 9 Be (dashed), d - 6 Li (full), 8 Li (dashed), e - Maxwellian 
fit for 7 Li, T s £ ope =8.9 MeV); (f) mean elemental event multiplicity N(Z) for QP fragments (solid 
point: experimental data, histogram: SMM predictions). 

The evidences of the equilibration of the QP source allow investigating some characteristics, 
such as the temperature and the excitation energy. We extract the temperature T through 
the method of double ratios of isotope yields [16]. As this method requires that the nuclei 
originate from the same emitting source, only the Z>3 fragments forward emitted (v par >6.5 
cm/ns, velocity of the QP in the laboratory frame) were considered. Table 1 gives the values 
for the most reliable thermometers [16,17]. Since experimental temperature measurements 
are affected by secondary decays, we also report in Table 1 the values corrected as suggested in 
Ref. [17]. Their mean value, T =3.9±0.3 MeV, can be considered as the break-up temperature 
of the QP decaying system. 



Table 1 

Table 1 

Temperatures of the QP system extracted from different double yield isotope ratio (T exp ) and 
calculated values after sequential feeding correction (T corr ). 



T exp (MeV) T corr (MeV) 



6 Li/ 7 Li - 


11 c/ i2 C 


3.9±0.2 


3.3±0.3 


9 Be/ w Be 


_11 C/ 12 C 


6.7±0.9 


3.7±0.5 


W B/ U B- 


.11 B /i2 B 


4.3±0.5 


4.3±0.5 


U B/ 12 B- 


- ll C/ 12 C 


4.3±0.3 


4.4±0.3 


U C/ 12 C- 


12 C/ 13 C 


3.7±0.2 


3.6±0.2 



Besides the indication given by the statistical model, we can give a rough estimate of the 
upper limit of the excitation energy using the energy conservation and assuming that on 
average there is an equal sharing of excitation energy between QP and QT. Then in the 
center of mass frame we have: Eq P = \{mpv 2 P — ttiqpVqp) + E* NECK , where nip, tuqp, 
vp, vqp are mass and velocity of projectile and QP, respectively. Neglecting the amount of 
energy transferred to the neck source (admittedly small to avoid its complete vaporization) 
and taking into account the mass difference between the projectile and the QP, the maximum 
of excitation energy ranges from 3.7 (if no nucleon transfer to a neck is assumed) up to 5.9 
MeV/nucleon (half nickel is lost in the reaction); if, for instance, we require the formation 
of an oxygen nucleus in the center of mass source we have accordingly an estimate of ~4.5 
MeV/nucleon. This rough estimation completely neglects other dissipation processes, as pre- 
equilibrium emission. 

The temperature and excitation energy values extracted in the present case are quite sim- 
ilar to those measured [2] for the QP fragmentation in the 197 Au + 197 Au reaction at 35 
MeV/nucleon (T =3.9±0.2 MeV, an upper limit of the excitation energy ~4.5 MeV/nucleon 
and a corresponding measured (and also predicted by SMM) value of ~4 MeV/nucleon). In 
Ref.[2] it has been found good agreement between the excitation energy values obtained 
following different methods (calorimetric [18], comparison with SMM predictions, rough es- 
timation of the upper limit), in the range (3-6 MeV/nucleon). 

In conclusion the present data analysis shows that the QP source has attained thermal 
equilibrium and that fragmentation is its main de-excitation process, well reproduced by a 
statistical approach. 

As already shown, in coincidence with the statistical fragmentation of the QP, we observed 
the emission of IMF from a IS. In the following we will discuss in detail the properties of 
this IMF emission at midrapidity. The expected distribution of the fragment velocity v par 
from the multifragmentation of the QP source is Gaussian. Thus the QP region of the v par 
distribution has been fitted, for each fragment Z value, with a Gaussian function, giving the 
experimental QP yield (Yqp). The fit parameters were extracted taking into account only 
the v par >6.5 cm/ns part. In Fig. 2 the results are plotted, together with the experimental 
yields. Similar behaviors have been found for all atomic numbers in the range 3<Z<14. Due 
to the experimental energy threshold (for the QT source side) the presented spectra are 
not symmetric around the cm. velocity. To avoid possible contamination in the midrapidity 
region we evaluated the yield (Yneck) from the IS source as twice the difference between 
the whole v par distribution for velocities higher than that of the cm. and Yq P . This was 



done to avoid distorsions due to efficiency effects and possible QT contaminations for the 
lowest velocities. We checked that in the considered Z range the experimental inefficiency 
doesn't affect the above-mentioned procedure [9]. 

In Fig. 4a the ratio between the relative yields (Y NE ck/Yqp) is presented as a function of 
the atomic number Z. We observe a bell-like shape, peaked around Z=9. We notice the high 
probability of IMF emission from the neck zone when the event IMF multiplicity is at least 
3 [4]. The presence of this maximum could be an effect of the breakup geometry; in fact, 
a rough consistency has been found with percolative calculations that compare emissions 
from a cylindrical shape neck, joining spherical QT and QP, and from the QP itself [5,14]. If 
on one side the QP disassembly is ruled by statistical models after thermal equilibrium has 
been reached, on the other side the neck emission exhibit quite different features that cannot 
be reproduced making statistical equilibrium assumptions. We have thus performed BNV 
calculations using different EOS parameters [3]. We found that with a compressibility term 
K of 200 MeV (soft EOS) there is an evident massive neck formation (after 200 fm/c), that 
is not reabsorbed by the QP or the QT (this behaviour disappears increasing the K values). 
These calculations predict that on average we have a Z=8 fragment in the neck zone, and 
show that the IS fragment production comes from material which is "surface-like" (since it 
originates from the overlap of the surfaces of the two nuclei) and which could be neutron rich 
[7]. To investigate the neutron content of the neck matter we observed the 6 He experimental 
yield. Its emission is quasi negligible from a statistical decay, both because the binding energy 
favours, for Z=2, a particle emission and because the N/Z value of 6 He is quite different from 
the corresponding ratio of the system (N/Z=2 for 6 ife,N/Z=30/28~l for the system). If the 
surfaces of the interacting nuclei are more neutron rich than the bulk matter or the dynamical 
process leading to the formation of a neck structure has a particular isospin dependence, the 
6 He production should be more abundant in the neck zone with respect to the QP zone. 
In Fig. 4b the experimental 4 He and 6 He yields are plotted versus v par '- one can see that, 
while the 4 He distribution is centered around the QP velocity (with a small emission from 
the neck zone) the 6 He yield is very scarse in the QP zone and starts to increase going 
towards the midvelocity region. One should take into account that the energy thresholds 
for mass identification are greater than for Z identification and that in the He case they 
produce cuts for velocities lower than 4 cm/ns. To better see the isospin effect in Figg.4c-d 
the Y( 6 He)/Y( A He) and Y( 3 He) /Y( A He) yield ratios are plotted. The amount of more rich 
(poor) neutron i/e-isotopes increases (decreases) going towards the midvelocity region. The 
behaviour of these ratios clearly shows the increase of neutron content in the midvelocity 
zone with respect to the interacting matter. This is a further indication that in midperipheral 
collisions we observe an IMF production which is due to two different mechanisms: one of 
statistical and the other of dynamical nature. 

In conclusion, in the study of the Ni+Ni 30 MeV/nucleon dissipative midperipheral collisions 
it has been possible to reveal events in which the IMF's are emitted by two different sources 
with different mechanisms. We are in presence of a QP (and a QT), with an excitation energy 
which leads to the multifragmentation regime; its decay can be fully explained in terms of 
a statistical disassembly of a thermalized system (T=3.9±0.2 MeV, E* ~4 MeV/nucleon). 
Contemporary to the IMF production from the QP source an intermediate source is formed, 
emitting both light particles and IMF. These fragments are more neutron rich than the 
average matter of the overall system and have a quite different charge distribution, with 
respect to the ones statistically emitted from the QP. These features can be considered as 
a signature of the dynamical origin of the midvelocity emission. The presented results then 
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Fig. 4. (a) Ratio of the measured yields for neck fragmentation and QP emission; (b) Yields 
of 4 He (dot-dashed line) and 6 He (full line, multiplied by a factor 10); Y( 6 He)/Y( A He) (c) and 
Y( 3 He)/Y( 4 He) (d) yield ratios as a function of Vpar- 

show that IMF can be produced via different mechanisms that can find contemporary room 
inside the same collision. 
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